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The current status of Ultra High Energy Cosmic Rays (UlfECR) is reviewed, with emphasis given to theoretical 
interpretation of the observed events. The galactic and extragalactic origin, in case of astrophysical sources of 
UHE particles, have the problems either with acceleration to the observed energies or with the fluxes and spectra. 
Topological defects can naturally produce particles with energies as observed and much higher, but in most cases 
fail to produce the observed fluxes. Cosmic necklaces and monopole-antimonopole pairs are identified as most 
plausible sources, which can provide the observed flux and spectrum. The relic superheavy particles are shown 
to be clustering in the Galactic halo, producing UHECR without Greisen-Zatsepin-Kuzmin cutoff. The Lightest 
Supersymmetric Particles are discussed as UHE carriers in the Universe. 
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1. Introduction 

Cosmic rays (CR) are observed in a wide en- 
ergy range, starting from subGeV energies and up 
to 3 • 10^0 eV (see Fig.l). Apart from the high- 
est energies, these particles are accelerated in our 
Galaxy, most probably, by shocks produced by 
SN II explosions. Up to energy lO^^ - lO^^ eV the 
CR flux is dominated by protons, at higher ener- 
gies CR have the mixed composition, and there 
are indications that at energies about ^ 10^^ eV 
iron nuclei dominate in the CR flux. In a wide 
range of energies from 1 GeV up the 3 • 10^^ eV 



the spectrum is power-law ^ E 



'2.65 



at energy 



3 • 10^^ eV the spectrum steepens and becomes 
- £'-31 at £ > 10" eV. At E > 10^^ eV a new 
more flat component appears (see Fig.l). The 
highest energies detected so far are 2 — 3-10^° eV. 

The first steepening of the spectrum (the knee) 
at energy 3 • 10^^ eV — 1 • 10^^ eV is usually ex- 
plained by inefficient confinement of CR in the 
Galaxy. This process must be accompanied by 
enrichment of heavy nuclei in CR composition at 
energy ~ 10^^ eV and above. There are some in- 
dications that such enrichment is really observed. 

There is no universal definition for Ultra High 
Energy Cosmic Rays (UHECR). Sometimes this 
term is apphed for £ > 1 • lO^*" eV or E > 1 • 
10^® eV. I shall use this term for £ > 1 • 10^^ eV, 
where the new flat component appears. 
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It is natural to think that this component has 
extragalactic origin, though, in principle, very 
large halo with regular magnetic field can con- 
fine particles of these energies, especially if they 
are heavy nuclei. UHECR of extragalactic ori- 
gin have a signature called the Greisen-Zatsepin- 
Kuzmin (GZK) cutoff 1^. This phenomenon is 
caused by energy losses of UHE protons due to 
pion production in collisions with microwave pho- 
tons. The energy losses start sharply increas- 
ing at £ ~ 3 • IQi^ eV (Fig.2). This energy is 
connected with energy of the spectrum steepen- 
ing ("cutoff') in the model-dependent way. In 
case the sources are distributed uniformly in the 
Universe (standard assumption), the steepening 
starts at Ebb ^ 3 ■ 10^^ eV. The flux at £ > Ebb 
is produced by nearby sources. If there is a local 
enhancement of the sources. Ebb increases ||];in 
case the sources are located at large distances, 
Ebb decreases and steepening is exponential. It 
is more convenient to characterize steepening by 
energy £1/2 where the flux becomes half of 
the power-law extrapolation of unmodified flux. 
In case of uniform distribution of the sources 
Ei/2 ~ 5.8 • 10"^^ eV [|j for a wide range of expo- 
nents 7 of generation spectrum. 

Apart from GZK cutoff, there may be two more 
signatures of extragalactic cosmic rays: a bump 
and a dip in differential spectrum which precede 
the cutoff The bump is a consequence of 

a number conservation of protons in the spec- 
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trum: protons loose energy and are accumulated 
before the cutoff. The dip is formed due to pair- 
production (e+e") energy losses of UHE proton. 
The both features show up most clearly in the 
differential spectrum of a single distant source in 
the case of a flat generation spectrum. In diffuse 
spectra (from many sources) these features are 
weak or absent. 

UHE nuclei spectra exhibit steepening (" cutoff' ) 
approximately at the same energy as protons, 
though due to different physical processes (see 
^ for a review). The relevant energy losses are 
caused by photodisintegration of nuclei at colli- 
sions with microwave photons, and the steepening 
energy is determined by energy, when photodis- 
integration energy-losses start to dominate over 
adiabatic ones (Fig. 2.). 

UHE photons with E-y - 10^^ - lO^^ eV have 
an absorption length less than 10 Mpc, mainly 
due to interaction with radio-background |^,^ . 

The observation of cosmic ray particles with 
energies higher than 10^" eV gives a serious chal- 
lenge to our understanding of origin of UHECR: 
What are the mechanisms of acceleration? Why 
the GZK cutoff is absent? 

2. Observational data 

The compilation ||^ of observational data for 
UHECR is given in Fig. 3. Two highest energy 
events correspond to energies 3 • 10^° eV (the 
Fly's Eye event §) and 2 • lO^o eV (the AGASA 
event ||l^). These energies are well above the 
GZK cutoff for uniformly distributed extragalac- 
tic sources. The particles with energies above 
10^° eV were observed in the past at the Hav- 
erah Park array |0| and Sydney array |12|. The 



Anisotropy of UHECR is not reliably observed. 
Some analyses (e.g. jl^ and jl^) indicate the 
excess of particles from Local Supercluster (LS) 
plane, while observations by AGASA [jl^ are con- 
sistent with a uniform distribution. On other 



latter detector has observed eight showers with 
E > 1 ■ lO^o eV. It operated from 1968 up to 
1979 and had large area, 87 km^. The scintilla- 
tor detectors were at 2m underground and hence 
only muon component of showers was measured. 
It is interesting to re-analyze the data using the 
new simulations for muon distribution |]l3|| . 
One shower with E 1.2 ■ lO^o eV was observed 
at the Yakutsk array ||l^,|l5|, though eight show- 
ers were expected if the Haverah Park data are 
correct. 



hand AGASA has observedy7| clustering of UHE 
events: three pairs of showers with angular sepa- 
ration less than 2.5° (for analysis see p^). 

Chemical composition as found from analysis 
of the Fly's Eye data |l^ is characterized by a 
change from predominantly heavy nuclei (iron) 
to the light nuclei at - 3 • 10^^ eV. The frac- 
tion of protons increases with energy and reaches 
90% at 10^^ eV The change of the chemical 
composition at £' 3 • 10^^ eV was not found in 
the muon data in AGASA experiment [Q. The 
data of Yakutsk array |Q also favors the proton 
composition at the highest energies. 



3. Galactic origin of UHECR 

There are two difficulties in an attempt to 
explain the observed events by sources in the 
Galaxy: the maximum acceleration energy is con- 
siderably less than 3 • 10^° eV and galactic mag- 
netic fields are too weak to isotropize UHE par- 
ticles. 

Acceleration to very high energies in the Galaxy 
can occur at the SN shocks, at galactic wind ter- 
minal shock and in young pulsars. 

The maximum energy for acceleration by 
shocks in the interstellar medium does not ex- 
ceed 10^^ eV It can be higher when the 
SN shock propagates through the region of pre- 
supernova stellar wind with strong magnetic field 
p3[ | , though the maximum energy is still less than 
3 • 10^^ eV Q. The galactic wind is expected 
to be terminated by a standing shock, where, 
in case of extreme values of parameters used in 
ref.(l^), particles can be accelerated up to ener- 
gies E ~ 10^" eV for iron nuclei. Most probably 
this energy is an order of magnitude less (see ^ 
and references therein). 

Another potential source of acceleration is 
young pulsar; in this case the maximum energy 
can in principle reach 10^^ eV However, in 
the concrete models of pulsar magnetosphere the 
maximum energy is less. 
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Propagation of UHECR in galactic magnetic 
fields was studied numerically in many works 
p9| . The crucial element of this analysis is pres- 
ence of regular magnetic field in galactic halo. All 
workers agree that at the highest observed ener- 
gies there must be very strong disc anisotropy, 
which obviously contradicts observations. 

We shall describe here the results of calcula- 
tions of ref.(||. 

Magnetic field is taken according to the model 
of ref.(|3^]). Several versions of the halo field is 
considered. The size of the halo is varied, but 
generically the large size of the halo from 10 to 30 
kpc is used. The sources of UHECR are assumed 
to be distributed uniformly in the disc. The tra- 
jectories of the antinuclei with rigidity R = E/Z, 
where Z is a charge of a nucleus, were followed 
step by step in the magnetic fields of the disc and 
halo. The flux in a given direction is proportional 
to the length of trajectory in the disc. An exam- 
ple of trajectories with rigidity 3 • 10^* V is shown 
in Fig. 4. The particles are emitted from the Earth 
in the different directions in Galactic plane z — 0, 
where z is the height over galactic plane. One can 
see that at energy i? ~ 3 • 10^* eV protons prop- 
agate almost rectilinear ly in the disc, producing 
thus strong disc anisotropy. In Fig. 5 the lifetimes 
of the particles in the disc, Td, and calculated 
anisotropy. A, are shown as a function of rigid- 
ity E/Z for two sizes of the halo. Since the flux 
is proportional to Td, one concludes from Fig. 5 
that at E > 3 ■ lO^^Z eV the particles typically 
do not return from the halo to disc. This con- 
clusion is confirmed by the lower part of Fig. 5, 
which shows large (disc) anisotropy at energies 
E > 3 ■ IQ^^Z eV. Note, that these calculations 
are performed for extreme case of very large mag- 
netic halo. 

Relativistic dust grains of galactic origin can 
be the carriers of UHE signal . The Lorentz 
factor of a dust grain should be large enough 
F > lO'^ to produce a nuclear-electromagnetic 
shower with muon component as observed in e.g. 
Sydney or AG AS A arrays. Approaching the sun 
such a grain accumulates the electric charge due 
to photoeffect produced by solar optical radiation 
(in the rest frame of a dust grain this is X-ray ra- 
diation). As a result a grain breaks up due to 



electric repulsion |^ . The dust grain hypothesis 
is also in contradiction with observed properties 
of UHE showers (J.Linsley). 

4. Extragalactic acceleration sources 

Shock acceleration (including ultra-relativistic 
shocks) and unipolar induction are the "stan- 
dard" acceleration mechanisms to UHE, consid- 
ered in the literature. These mechanisms can op- 
erate in the various astrophysical objects, such as 
Active Galactic Nuclei (AGN), large scale struc- 
tures (e.g. the shocks in AGN jets or shocks in 
the clusters of galaxies), in gamma-ray bursters 
(ultra-relativistic shocks), in the accretion discs 
around massive black holes (due to large elec- 
tric potentials produced by unipolar induction) 
etc. A comprehensive list of possible sources was 



recently thoroughly studied in ref.(|33|) with a 
conclusion, that maximum energy of acceleration 
does not exceed 10^^ — 10^° eV. The most promis- 
ing source from this list is a hot spot in radio- 
galaxy produced by a jet p^-|3^. 

A powerful jet ejected from the AGN supplies 
energy to a gigantic radiolobe. A hot spot ob- 
served at the termination of the jet is interpreted 
as a location of a standing shock. This is an ideal 
place for acceleration of protons to very high en- 
ergies: magnetic field is strong and the energy 
density of radiation, responsible for proton energy 
losses, is relatively small. The maximum energy 
can be estimated as [B4| 



E„ 



1 • W"H^4{R/1 kpc)vj eV, 



(1) 



where H^^ is the magnetic field in units of 
10"'* G, R is radius of the shock and Vj is velocity 
of the jet in units of sound speed. However, the 
powerful radio sources are at large distances from 
our Galaxy, and the maximum energy is strongly 
attenuated. The discussed sources can provide 
the observed flux up to energies 6 — 7 • 10*^ eV. 

Unipolar induction produces very large poten- 
tial drop in the accretion discs around massive 
black holes (see ref.(^,Q and also |^ for a dis- 
cussion). The electrical potential for a rotating 
disc, at the distance r from a black hole is 



(r) = -^HcrJnR/r 
V6 



(2) 
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where Tc = 9 • lO'^iMh/MQ) is the radius of the 
last stable orbit for a black hole with mass Mh, 
He is magnetic field at the last stable orbit, and 
i? is a radius of accretion disc. The maximum 
potential given by Eq.d) is 0c - 3 • 10^1 V for 
Mh ^ 1 • lO^M© and ff c ~ 1 • 10'' G. This mecha- 
nism is attractive because it can operate not only 
in AGN, but also for the old black holes, which 
lost their activity (M in the disc is small). Such 
sources can be located nearby, e.g. in the Local 
Supercluster, from where UHE protons can reach 
us without appreciable energy losses. 

Relativistic shocks can in principle provide very 
high maximum energy . A particle reflecting from 
a relativistically moving mirror increases its en- 
ergy by a factor proportional to F^, where T is 
a Lorentz factor of the mirror. This acceleration 
phenomenon is known from the time of the pi- 
oneering work by E.Fermi p9[ |. The less known 
phenomenon is capturing of accelerated particles 
behind relativistically moving shock front. One 
can easily reconstruct this principle considering 
head-on collision of a particle with a transverse 
relativistic shock (magnetic field is perpendicu- 
lar to the shock normal). Let us assume that 
magnetic field behind the shock is homogeneous 
on the scale of the particle Larmor radius. Then 
at the moment a particle finishes semi-circle, the 
shock front run away to a distance cTl/2 from 
a particle, where Tl is the Larmor period. The 
confinement described above can be illustrated by 
numerical simulation in ref.([Q) shown in Fig. 6. 
In Fig. 6a the shock velocity is non-relativistic 
(0.1c) and a particle refiects many times from a 
shock. In Fig. 6b the shock is relativistic, 0.87c, 
and a particle is captured behind the shock. The 
same phenomenon is clearly seen in the relativis- 
tic transverse shocks in electron-positron plasma, 
( Fig. 4 from |41|). Somewhat more compli- 



cated, but similar mechanism operates for par- 
allel shock, i.e. when magnetic field is parallel to 
the shock velocity. 

The capturing mechanism described above, 
does not exclude completely the F^-regime of ac- 
celeration at relativistic shocks; it restricts the in- 
cident angles at which particles escape and thus 
the flux of accelerated particles. The most inter- 
esting objects where F^-mechanism might operate 



are gamma-ray bursters 43|. The Lorentz fac- 
tor of the shock here can reach 10^ - 10^. How- 
ever, the capturing properly taken into account 
might dramatically decrease the output of accel- 
erated particles. As to the explanation of the ob- 
served UHECR, the protons from cosmologically 
remote gamma-ray bursts strongly degrade in en- 
ergy on the way to our Galaxy. 

Astrophysical sources of observed UHECR 
must satisfy the observational constraints. The 
absence of GZK cutoff at i?i/2 ~ 6 • 10^^ eV con- 
tradicts the uniform distribution of the sources in 
extragalactic space. If the sources are located as 
a dense group around our Galaxy, -E1/2 increases 
with increasing of density contrast, i.e. the ratio 
of number density of sources inside the group and 
outside it [|j . The Local Supercluster LS) can re- 
alize this possibility. This model was developed 
in ref.(@). The UHECR sources in LS can be 
the old massive black holes with large electric po- 
tential induced in the accretion discs by unipolar 
induction. For increasing E1/2 up to 1 • 10^" eV 
the density contrast larger than 10 is needed. The 
anisotropy can be smaller than the observed one 
if magnetic fleld in superclusters is as large as 
10"^ G 11. 

5. Topological defects and relic particles. 

Topological defects, TD, (for a review see p^]) 
can naturally produce particles of ultrahigh ener- 
gies (UHE). The pioneering observation of this 
possibility was made by Hill, Schramm and 



Walker |47| (for a general analysis of TD as UHE 
CR sources see |Q and for a recent review [p9|). 

In many cases TD become unstable and de- 
compose to constituent flelds, superheavy gauge 
and Higgs bosons (X-particles), which then de- 
cay producing UHECR. It could happen, for ex- 
ample, when two segments of ordinary string, or 
monopole and antimonopole touch each other, 
when electrical current in superconducting string 
reaches the critical value and in some other cases. 

In most cases the problem with UHECR from 
TD is not the maximal energy, but the fluxes. 
One very general reason for the low fluxes con- 
sists in the large distance between TD. A dimen- 
sion scale for this distance is the Hubble distance 
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Hq^. However, in some rather exceptional cases 
this dimensional scale is multiplied to a small di- 
mensionless value r. If a distance between TD 
is larger than UHE proton attenuation length, 
then the flux at UHE is typically exponential sup- 
pressed. 

Ordinary cosmic strings can produce particles 
when a loop annihilate into double line ||5^. The 
produced UHE CR flux is strongly reduced due 
to the fact that a loop oscillates, and in the pro- 
cess of a collapse the two incoming parts of a 
loop touch each other in one point producing thus 
the smaller loops, instead of two- line annihilation. 
However, this idea was recently revived due to re- 
cent work . It is argued there that the energy 
loss of the long strings is dominated by produc- 
tion of very small loops with the size down to the 
width of a string, which immediately annihilate 
into superheavy particles. A problem with this 
scenario is too large distance between strings (of 
order of the Hubble distance). For a distance be- 
tween an observer and a string being the same, 
the observed spectrum of UHE CR has an expo- 
nential cutoff at energy E ^ 3 ■ 10^^ eV. 

Superheavy particles can be also produced 
when two segments of string come into close con- 
tact, as in cusp events [Q. This process was 
studied later in ref.(|5^) with a conclusion that 
the resulting cosmic ray flux is far too small. An 
interesting possibility suggested in ref.(|^^) is the 
cusp "evaporation" on cosmic strings. When the 
distance between two segments of the cusp be- 
comes of the order of the string width, the cusp 
may "annihilate" turning into high energy parti- 
cles. , which are boosted by a very large Lorentz 
factor of the cusp However, the resulting 

UHE CR flux is considerably smaller than one 
observed [ p4[ . 

Superconducting strings p5[ | appear to be much 
better suited for particle production. Moving 
through cosmic magnetic fields, such strings de- 
velop electric currents and copiously produce 
charged heavy particles when the current reaches 
certain critical value. The CR flux produced 
by superconducting strings is affected by some 
model-dependent string parameters and by the 
history and spatial distribution of cosmic mag- 
netic fields. Models considered so far failed to 



account for the observed flux . 

Monopole-antimonopole pairs (MM) can form 
bound states and eventually annihilate into UHE 
particles |^^, js^. For an appropriate choice of 
the monopole density um, this model is consistent 
with observations; however, the required (low) 
value of ua/ implies fine-tuning. In the first phase 
transition G ^ H x U{1) in the early Universe 
the monopoles are produced with too high den- 
sity. It must then be diluted by inflation to very 
low density, precisely tuned to the observed UHE 
CR flux. 

Monopole- string network can be formed in 
the early Universe in the sequence of symmetry 
breaking 



H X U{1) H X 



N- 



(3) 



For > 3 an infinite network of monopoles 
connected by strings is formed. The magnetic 
fluxes of monopoles in the network are channeled 
into into the strings that connect them. The 
monopoles typically have additional unconfined 
magnetic and chromo-magnetic charges. When 
strings shrink the monopoles are pulled by them 
and are accelerated. The accelerated monopoles 
produce extremely high energy gluons, which 
then fragment into UHE hadrons ||5^]. The pro- 
duced flux is too small to explain UHE CR ob- 
servation Jsot . 

Cosmic necklaces are TD which are formed in 
a sequence of symmetry breaking given by Eq. (|^) 
when N = 2. The first phase transition produces 
monopoles, and at the second phase transition 
each monopole gets attached to two strings, with 
its magnetic flux channeled along the strings. The 
resulting necklaces resemble "ordinary" cosmic 
strings with monopoles playing the role of beads. 
Necklaces can evolve in such way that a distance 
between monopoles diminishes and in the end all 
monopoles annihilate with the neighboring anti- 
monopoles Q. 

An important quantity for the necklace evolu- 
tion is is the dimensionless ratio r = m/ ^d, where 
m is the monopole mass, fi is is the string ten- 
sion, determined by U{1) symmetry scaling scale 
rjs , and d is a distance between two neighbouring 
monopoles. In ref.(|61|) it is argued that in the 
process of the necklace evolution r(t) is driven 
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towards large values r 3> 1. The characteristic 
length scale, equal to the typical separation of 
necklaces, 

e - r~^^% (4) 

is much smaller than horizon t, when r becomes 
large. 

A requirement for all models explaining the ob- 
served UHE events is that the distance between 
sources must be smaller than the attenuation 
length for UHE particles. Otherwise the flux at 
the corresponding energy would be exponentially 
suppressed. This imposes a severe constraint on 
the possible sources. For example, in the case of 
protons with energy E ^ (2 — 3) • 10^" eV, the 
proton attenuation length is 19 Mpc . If protons 
propagate rectilinearly, there should be several 
sources inside this radius; otherwise all particles 
would arrive from the same direction. If particles 
are strongly deflected in extragalactic magnetic 
fields, the distance to the source should be even 
smaller. Therefore, the sources of the observed 
events at the highest energy must be at a dis- 
tance i? < 15 Mpc in the case or protons. 

For the necklaces the distance between sources, 
given by Eq.(^, satisfies this condition for r > 
3 • 10"*. This is in contrast to other potential 
sources, including supeconducting cosmic strings 
and powerful astronomical sources such as AGN, 
for which this condition imposes severe restric- 
tions. 

The diffuse fluxes of UHE protons and photons 
from necklaces are given in Fig. 7. One can see 
that at E' ~ 1 • 10^" eV the gamma-ray flux is 
considerably lower than that of protons. This 
is mainly due to the difference in the attenua- 
tion lengths for protons (110 Mpc) and photons 
(2.6 Mpc and 2.2 Mpc §1). At higher energy 
the attenuation length for protons dramatically 
decreases (13.4 Mpc at E = 1 ■ 10^^ GeV) and 
the fluxes of protons and photons become com- 
parable [|6l| . 

The predictions in Fig. 7. are compared with 
the AGASA data |^^. The agreement is rather 
good for E > 2 ■ 10^^ eV. The contribution of 
UHE photons increases the total flux at > 
2 • 10^° eV. The contribution of low-energy com- 
ponent (e.g. from radiogalaxies [|^) can easily 



improve the agreement at lower energies. 

The radiation from Topological Defects can ex- 
plain the diffuse gamma-ray background above 
10 GeV §|. 

Superheavy relic particles can be sources of 
UHE CR [|4l-|§. In this scenario Cold Dark 
Matter (CDM) have a small admixture of long- 
lived superheavy particles. These particles must 
be heavy, mx > 10^^ GeV, long-lived tx > to, 
where to is the age of the Universe, and weakly 
interacting. The required life-time can be pro- 
vided if this particle has (almost) conserved quan- 
tum number broken very weakly due to warmhole 
p5| | or instanton ||6^ effects. Several mechanisms 
for production of such particles in the early Uni- 
verse were identified. Like other forms of non- 
dissipative CDM , X-particles must accumulate 
in the halo of our Galaxy and thus they pro- 
duce UHE CR without GZK cutoff and without 
appreciable anisotropy. 

The more detailed discussion as well as cal- 
culated fluxes of UHE protons and photons are 
given in the paper by M. Kachelriess (these Pro- 
ceedings) . 

The characteristic and unavoidable feature of 
this model is an excess of gamma-ray flux over the 
nucleon flux. It follows from the more effective 
production of pions than nucleons in the QCD 
cascades from the decay of X -particles and from 
absence of absorption in the Galactic halo. 

The spectrum of the observed FAS is formed 
due to fluxes of gamma-rays and nucleons. The 
gamma-ray contribution to this spectrum is 
rather complicated. In contrast to low energies, 
the photon-induced showers at i? > 10^* eV 
have the low-energy muon component as abun- 
dant as that for nucleon-induced showers ]6^ . 
The shower production by the photons is, in prin- 
ciple, suppressed by the LPM effect and by 
absorption in geomagnetic fleld. However, as was 
noted in cascading in geomagnetic field re- 
sults in arrival of a bunch of photons (each with a 
smaller energy than the primary one) at the top 
of atmosphere. They produce one shower with 
LPM effect reduced because of the smaller ener- 
gies of photons in the bunch (see |6^] and refer- 
ences therein for further discussion). Experimen- 
tal discrimination of gamma-ray induced showers 
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from that produced by protons is very important 
task. 

The UHE carriers can be, in principle, not only 
protons and photons, but other particles, such 
as neutrinos, gluinos [[70|-[72|, neutralinos |^2| and 
monopoles |7^j7l[]. In the next section we shall 
shortly discuss discuss the Lightest Supersym- 
metric Particle as a candidate for a carrier of a 
signal in UHECR detectors. 

6. LSP as UHE carrier 

The Lightest Supersymmetric Particle (LSP) 
can be either stable, if R-parity is strictly con- 
served, or unstable, if R-parity is violated. To 
be able to reach the Earth from most remote re- 
gions in the Universe, the LSP must have lifetime 
longer than tlsp > ^o/L, where to is the age of 
the Universe and F = E /vtllsp is the Lorentz- 
factor of the LSP. In case niLSP ^ 100 GeV, 
TLSP > 1 yr- 

Theoretically the best motivated candidates for 
LSP are the neutralino and gravitino; the latter 
is practically undetectable as UHE particle. 

In all elaborated SUSY models the gluino is 
not the LSP. From experimental point of view 
there is some controversy if the low-mass win- 
dow 1 GeV < rrig < 4 GeV for the gluino is 



still allowed [ [74[|75| ] . Recently the light-gluino was 
claimed to be ruled out on the basis of its con- 
tribution to (3 function for the strong interaction 
76 1 . Finally, there is one more argument about 
light gluino ||7^ : If gluino is stable or quasi-stable 
and if the lightest gluino-baryonic state is guud 
, this heavy hydrogen is overproduced by cosmic 
rays in the Earth atmosphere. Nevertheless, we 
shall discuss the gluino as UHE carrier |f70|-[72[. 
We shall refer to any colorless hadron contain- 
ing gluino as ^-hadron. The lightest ^-hadron 
is most probably glueballino gg. It is stable if 
gluino is LSP. The lightest gluino-baryonic state, 
gluebarino, is almost stable because of very weak 
violation of baryonic number. Light glueballinos 
as UHE particles with energy E > 10^^ eV were 
considered in some detail in the literature in con- 
nection with Cyg X-3 

UHE LSP are most naturally produced at the 
decays of unstable superheavy particles, either 



from TD or the relic ones ||7^ . 

The QCD parton cascade is not a unique cas- 
cade process. A cascade multiplication of partons 
at the decay of superheavy particle appears when- 
ever a probability of production of extra parton 
has the terms a In or a In^ , where Q is a 
maximum of parton transverse momentum. Re- 
gardless of smallness of a, the cascade develops as 
far as a In (5^ > 1. Therefore, for extremely large 
we are interested in, a cascade develops due to 
parton multiplication through SU{2) x C/(l) in- 
teractions as well (see |72 for calculations). LSP 
take away a considerable fraction of the total en- 
ergy (~ 40%). 

Neutralino as UHE carrier has no much advan- 
tages over neutrino: the neutralino fluxes, pro- 
duced at the decay of superheavy particles, are 
less than neutrino fluxes and neutralino-nucleon 
cross-section at very high energy is less than that 
for neutrino [ [72| . 

Light gluino is effective as the UHE carrier. 
The energy losses of glueballino on a way from a 
source to the Earth is less than for a proton. The 
dominant energy loss is due to pion production in 
collisions with microwave photons. Pion produc- 
tion effectively starts at the same Lorentz-factor 
as in the case of the proton. This implies that 
the energy of the GZK cutoff is a factor rrig/mp 
higher than in case of the proton. The attenu- 
ation length also increases because the fraction 
of energy lost near the threshold is small, ji/mg, 
where /i is a pion mass. Therefore, even for light 
glueballino, rrig > 2 GeV , the steepening of the 
spectrum is less pronounced than for protons [ f72[ . 

A very light UHE glueballino interacts with the 
nucleons in the atmosphere similarly to UHE pro- 
ton. The cross-section for heavy glueballino with 
nig > 150 GeV is small for large energy transfer 
needed for production of extensive air showers, 
and thus these particles cannot be responsible for 
the observed UHE events. 

Thus, only UHE gluino from the low-mass win- 
dow 1 GeV > rUg < 4 GeV could be a candidate 
for observed UHE particles, but it is disfavored 
by the arguments given above. 
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7. Conclusions 

At £^ > 1 • 10^^ eV a new component of cos- 
mic rays with a flat spectrum is observed. Two 
highest energy events have i? w 2 — 3 • 10^" eV. 
According to the Fly's Eye and Yakutsk data the 
chemical composition is better described by pro- 
tons than heavy nuclei. The AGASA data are 
consistent with isotropy in arrival of the particles, 
though theoretical analysis reveals some correla- 
tion of arrival direction with Local Supercluster 
plane. AGASA has observed clustering of UHE 
events: three pairs of particles with small angular 
separation. 

The galactic origin of UHECR is disfavored: 
the maximal observed energies are less than that 
known for the galactic sources, and the strong 
Galactic disc anisotropy is predicted even for the 
extreme magnetic fields in the disc and halo. 

The signature of extragalactic UHECR is GZK 
cutoff. The position of steepening is model- 
dependent value. For the Universe uniformly 
filled with sources, the steepening starts at Ehb ~ 
3 • IQis eV and has « 6 • lO^^ eV (the energy 
at which spectrum becomes a factor of two lower 
than a power-law extrapolation from lower ener- 
gies). The spectra of UHE nuclei exhibit steepen- 
ing approximately at the same energy as protons. 
UHE photons have small absorption length due 
to interaction with radio background radiation. 

The extragalactic astrophysical sources theo- 
retically studied so far, have either too small 
Emax or are located too far away. The Local 
Supercluster (LS) model can give spectrum with 
Ei/2 ^ 10^° eV, if density contrast for the sources 
(the ratio of densities inside LS and outside) is 
larger than 10. 

Topological Defects naturally produce particles 
with extremely high energies, much in excess of 
what is presently observed. However, the fluxes 
from most known TD are too small. So far only 
necklaces and monopole-antimonopole pairs can 
provide the observed flux of UHE CR. 

Another promising sources of UHE CR are relic 
superheavy particles. These particles should be 
clustering in the halo of our Galaxy, and thus 
UHECR produced at their decays do not have 
the GZK cutoff. The signatures of this model are 



dominance of photons in the primary flux and 
Virgo cluster as a possible discrete source. 

Apart from protons and photons, the light 
gluinos can be successful UHE carriers, but they 
are disfavored in mass interval at interest. 
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Figure 1: Compilation of cosmic ray spectrum from [1]. Spectrum of UHECR is taken according 
to AGASA data only, because the spectrum shape is better seen from the data of one array. 
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gure 2: Energy losses of UHE protons and nuclei [3] 
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Figure 3: Compilation of cosmic ray spectrum from [1]. Spectrum of UHECR is taken according 
to AGASA data only, because the spectrum shape is better seen from the data of one array. 
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Figure 4: Trajectories of antinuclei with rigidity E/Z = 3 • 10^^ eV emitted in the different 
directions in the Galactic plane {z = 0). The upper figure gives projection of the trajectories 
on the galactic plane {x,y); two lower figures - projections on the planes {z,x) and {z,y), 
perpendicular to the Galactic plane. The distances are in kpc. 
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This figure "figS.jpeg" is available in "jpeg" format from: 



http://arXiv.org/ps/hep-ph/9802351vl 




Figure 6: Capturing of accelerated particles behind the front of a relativistic shock [41]. Fig- 
ure 6a corresponds to non-relativistic shock, Vgh = 0.1c, and Fig. 6b to relativistic shock 
Vsh = 0.87c. The position of the shock is shown by a dashed line, the direction of the ve- 
locity by an arrow. 
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Figure 7: UHE proton flux from cosmic necklaces with r^^u ^ 1 • 10^^ GeV^ and 
nix ^ 1 • 10^'' GeV. Curve with p shows the proton flux, 7 is for photon flux [62]. 
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